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ABSTRACT 

In this paper we examine how well galaxies and intra-cluster gas trace the gravita- 
tional potential of clusters. Utilizing mass profiles derived from gravitational lensing 
and X-ray observations, coupled with measured galaxy velocities, we solve for the 
velocity anisotropy parameter /3 or b(r) using the anisotropic Jeans equation. This is 
done for five clusters, three at low redshift: A2199, A496 and A576 and two at high 
redshifts: A2390 and MS1358. With X-ray temperature profiles obtained from Chan- 
dra and ASCA/ROSAT data, we estimate flx(r) the ratio of energy in the galaxies 
compared to the X-ray gas. We find that none of these clusters is strictly in hydro- 
static equilibrium. We compare the properties of our sample with clusters that form in 
high-resolution cosmological N-body simulations that include baryonic physics. Sim- 
ulations and data show considerable scatter in their (3 OT b(r) and (3x{r) profiles. We 
demonstrate the future feasibility and potential for directly comparing the orbital 
structure of clusters inferred from multi-wavelength observations with high resolution 
simulated clusters. 

Key words: 



1 INTRODUCTION 

Galaxy dynamics in clusters provided the first evidence of 
dark matter (Zwicky 1933; 1937). However, the orbital dis- 
tribution of galaxies remains a major source of uncertainty 
in deriving mass profiles from cluster galaxy dynamics. Sev- 
eral authors have shown that a wide variety of mass profiles 
can reproduce the observed galaxy and velocity distribution 
in clusters (The & White 1986; Merritt 1987). The obser- 
vational errors in the measurements of velocity dispersion 
profiles can translate to large uncertainties in dynamical 
mass estimates. Even when the kinematic data available are 
extensive, simplifying assumptions are usually necessary to 
estimate orbital distributions from velocity data alone (van 
der Marel 2000; Biviano & Katgert 2004). The key assump- 
tion in using galaxy kinematics to derive the gravitational 
potential of the dark matter is that they are robust tracers 
of the total mass distribution. However, as galaxies repre- 
sent only a small percentage (~ 1%) of the total mass in 
clusters, their velocity distribution could in principle differ 
significantly from that of the dark matter. Here we explore 
the feasibility of using independent mass estimators to break 



degeneracies in dynamical models. We can then measure the 
orbital distribution of cluster galaxies without assuming that 
cluster galaxies are good tracers of the velocity field. 

If as expected in the standard scenario, galaxies form 
primarily at the centers of dark matter halos, the orbital dis- 
tribution of cluster galaxies is likely to be closely connected 
to that of the dark matter. However, various physical mech- 
anisms could produce a velocity bias between the two com- 
ponents. For instance, the formation of galaxy pairs could 
transform orbital energy into internal energy and thus pro- 
duce velocity anti-bias (Fusco-Femiano & Menci 1995). Dif- 
ferent investigations with numerical simulations have pro- 
duced a range of predictions for the magnitude and even 
the sense of velocity bias (Colin, Klypin & Kravtsov 2000, 
Ghigna et al. 2001, Diemand, Moore & Stadel 2005, Gao 
et al. 2005, Faltenbacher et al. 2005). Colm et al. (2000) 
find that cluster galaxies are positively biased in the ra- 
dial range (0.2-0.8)r v i r , where r v i r is the virial radius, and 
approximately unbiased at larger radii. They attribute this 
radial dependence to dynamical friction on cluster galaxies 
near the cluster centre: dynamical friction is more efficient at 
lower velocities, so these galaxies merge, and the remaining 
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galaxies are therefore preferentially high-velocity galaxies. 
Faltenbacher et al. (2005) find that galaxies move slightly 
faster than the dark matter particles with a velocity bias 
factor of approximately 1.1. 

Gravitational lensing offers perhaps the cleanest and 
most powerful way to map mass distributions of galaxy clus- 
ters. The lensing mass estimates can be used directly to 
solve for the orbits of galaxies in the cluster. This approach 
was proposed and outlined in Natarajan & Kneib (1996) for 
clusters with mass estimates derived by combining strong 
and weak lensing effects. Their proposed formalism can be 
applied to any estimator of the mass profile that is inde- 
pendent of the velocity structure and therefore be used to 
robustly constrain the velocity anisotropy parameter. One 
such estimator is X-ray observations of the temperature and 
surface brightness distribution of the intra-cluster medium 
(ICM). In this paper, using X-ray and lensing mass profiles 
we study five clusters divided into two samples, three clus- 
ters (A2199, A496 and A576) from the Cluster And Infall 
Region Nearby Survey (CAIRNS; Rines et al. 2003) that 
constitute the nearby sample and two clusters (A2390 and 
MS1358) that form our distant sample 1 . We solve for the 
nature of galaxy orbits via the velocity anisotropy parame- 
ter (3 OI b{r), using the measured galaxy velocities, the line-of- 
sight velocity dispersions and radial mass profiles. Recently, 
Mathews & Brighenti (2003) in a similar vein used the X- 
ray temperature profile of the elliptical galaxy NGC 4472 to 
constrain the orbital distribution of its stars. 

There are likely to be systematic uncertainties in both 
X-ray and lensing mass estimates. For instance, the lens- 
ing signal around a cluster may have contributions from 
associated large-scale structure, and the overall normaliza- 
tion is typically unknown due to mass-sheet degeneracy. The 
mass-sheet degeneracy can however be broken by the com- 
bined use of strong and weak lensing features to constrain 
the mass profile (Natarajan et al. 2002; Bradac, Lombardi 
& Schneider 2004). On the other hand, the intracluster gas 
may not be in equilibrium; Chandra observations show evi- 
dence of buoyant bubbles near the cores of clusters (McNa- 
mara et al. 2000; Fabian et al. 2000; Churazov et al. 2001) as 
well as cold fronts and substructure (Markevitch et al. 2000; 
Vikhlinin, Markevitch, & Murray 2001). However, applying 
these techniques to clusters has an advantage over ellipticals 
because the hot gas is the dominant component of baryonic 
mass in clusters, whereas it is a much smaller component 
of the baryonic mass in ellipticals. For the clusters in our 
sample, the lensing mass estimates have been obtained us- 
ing strong and weak lensing data enabling the calibration 
of the mass, therefore the mass-sheet degeneracy has been 
lifted (Natarajan et al. 2005). 

The combined X-ray and optical data can constrain the 
properties of both baryonic tracers - the galaxies and the 
gas. One of these constraints comes from the so-called fix 
problem 2 . The relative internal energies of galaxies and gas 

1 Note that we use the word sample to define our chosen clusters, 
which span a range of morphology, mass and X-ray luminosity. 

2 Unfortunately, the traditional notation uses /3 as the symbol for 
both the orbital anisotropy of tracers and the ratio of energy per 
unit mass in galaxies and gas. We use subscripts to distinguish 
the two uses: /3 or b refers to orbital anisotropy and fix refers to 
the energy ratio. 



in clusters can be characterized by the parameter fix = 
prnipdl I (kTx) where fi is the mean molecular weight, m p is 
the proton mass, k is Boltzmann's constant, Tx is the X-ray 
temperature of the ICM, and <r 2 is the radial velocity disper- 
sion (usually estimated by assuming isotropic orbits so that 
<t 2 equals the line-of-sight velocity dispersion). The relative 
spatial distributions of galaxies and gas can be characterized 
by a different parameter, usually called f3 ga i, which should 
equal fix if the populations are in equilibrium. 

The ICM is assumed to be a fairly good tracer of the 
cluster potential, despite the fact that it is unlikely to be in 
strict hydro-static equilibrium. We are now able to test the 
robustness of this assumption, as we can compare the mea- 
sured temperature profile to that predicted from the orbital 
model for the cluster. 

There is growing evidence that the ICM is not in ther- 
mal equilibrium, in that instance it is unclear how to inter- 
pret fix- Many studies find an average value of (5x ~ 0.67 
(Jones & Forman 1984, 1999). Values of fix < 1 suggest 
that the ICM contains more energy per unit mass than the 
galaxies. This situation could arise from energy input in the 
ICM or by some of the orbital energy of the galaxies being 
converted into internal energies (Fusco-Femiano & Menci 
1995). We use either X-ray or lensing mass profiles to de- 
termine the orbital distribution /3 or b(r), which then allows 
us to estimate a r (r) and hence f3x(r), the energy ratio as 
a function of radius. Additionally, we compare the results 
for our clusters with high-resolution cosmological N-body 
simulations that include hydrodynamics (Faltenbacher et al. 
2005). These simulations trace the evolution of the gas and 
dark matter and offer the ideal comparison set for our di- 
verse clusters. 

In this paper we demonstrate that with the wealth of 
multi- wavelength data likely to be available in the very near 
future for large cluster samples, fruitful comparisons can be 
made with samples of simulated clusters selected in a simi- 
lar fashion from high-resolution cosmological N-body simu- 
lations. All the observational data to do so are available at 
the present time only for a handful of clusters that we study 
here. 

The outline of our paper is as follows: in Section 2, we 
present the formalism and the framework for our analysis. 
In Section 3, we discuss the derivation of the mass profiles 
from lensing and X-ray observations. We discuss the galaxy 
velocity data in Section 4. We present the derived orbital 
distributions in Section 5. Section 6 discusses the entropy 
profiles of the ICM. We compare our calculations with sim- 
ulations in Section 7 and we conclude with a discussion of 
the results in Section 8. We assume a cosmology of Q. m = 0.3, 
Qa = 0.7, and H = 70 km s _1 Mpc -1 throughout. 

2 THE ANISOTROPIC JEANS EQUATION 

The velocity anisotropy profiles of galaxy clusters can be 
obtained by solving the anisotropic Jeans equation, 

d(z/ g a 2 ) 2/3 orb (r)^ g o- 2 _ GM to t(r)v g , . 

dr r r z 

which applies to collisionless, spherically symmetric systems 
of particles. Here M tot (r) is the total cluster mass contained 
within radius r as determined by gravitational lensing or 
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X-ray observations, v g {r) is the three-dimensional galaxy 
number density, crf.{r) is the radial velocity dispersion of 
the galaxies, and f3 or b(r) is the velocity anisotropy . This 

CT 2 

parameter /? rb(r) is defined to be (1 — and is a measure 
of the predominance of tangential orbits over radial ones. 
For isotropic orbits /3 or b = 0; < /3 or b < 1 for radial orbits 
and for tangential orbits — oo < /3 OI b < 0. The anisotropic 
Jeans equation above can be combined with the equation 
that defines the measured line-of-sight velocity dispersion 
aios(R), 



±[Z g (R)aL(R)} = 



- R 



ru g (r)a 2 (r)dr 
R Vr 2 - R 2 ~ 

f3 olb (r)a 2 (r)p g (r)dr 



f 

Jr. 



rVr 2 - R 2 



(2) 



to yield two independent integro-differential equations for 
<r 2 (r) and (3 or b(r). In the latter equation, S 9 (i?) is simply 
the projected two-dimensional galaxy number density, which 
can be obtained by integrating v g (r). Following the proce- 
dure outlined in Natarajan & Kneib 1996 (and in Bicknell 
et al. 1989), one can solve these equations to express a 2 (r) 
as a sum of four integrals, 

v g (r)ol = Tcrmi(r) - Term 2 (r) + Term 3 (r) - Term 4 (r), (3) 

f-Rt 



Termi(r) = - 



GMtotir)^, 

-9. U ' ' 



Term 2 (r) = — / GM to t{r)v g rdr, 



Term3 ( 



RE g (R)o los (R)dR, 



(4) 
(5) 
(6) 



cflt 



2 r 

Term 4 (r) = — / RE g (R)af os (R) 

J r 

( - — sin -1 — ) dR. 

V JB 2 - r 2 R) 



(7) 



Here Rt is a large radius at which both f g (r) and S 3 (r) 
asymtote to zero. We find that using R t = 3.5 Mpc re- 
sults in smooth /3 or b (r) profiles in the range — 1 Mpc. The 
above integrals can easily be computed numerically using a 
constant step size trapezoid-rule algorithm. In the case of 
Terrru, a simple change of variable, R — r-v/l + C 2 > is used 
to eliminate the integrable singularity at the lower limit. 
After computing a 2 , one can obtain /3 rb(r) from the Jeans 
equation, 



>(»•) = -; 



GMtot{r)v g 



+ —{v g a r ) 



(8) 



3 CLUSTER MASS PROFILES 

To solve the anisotropic Jeans equation above, we use the in- 
dependently determined mass profile for clusters M{r) from 
lensing and/or X-ray data. The total mass profiles of all 

3 The reader should note that all quantities used in this paper are 
three-dimensional, unprojected quantities, unless explicitly stated 
otherwise. 



the clusters are fitted with the NFW model. This 'universal' 
profile provides the best-fit to the structure of dark mat- 
ter halos that form in N-body simulations in the context 
of a cold dark matter cosmogony (Navarro, Frenk & White 
1996). The NFW density profile can be written as 



P(r) = $c 

pcrit (r/r s )(l + r/r s ) 2 



(9) 



where p(r) is the three-dimensional total density, p cr it = 
3H 2 /8tvG is the critical density (H is the value of the Hubble 
constant at the cluster's redshift), r s is the characteristic 
scale radius, and 



200 



3 pn(l + c) - c/(l + c)] ' 



(10) 



The concentration parameter is defined by c = r2oo/r 3 , 
where r2oo is the radius within which the average density is 
200 times the critical density. The normalization, and hence 
the concentration c, of the NFW model fits for our distant 
clusters is derived from the detected strong lensing features 
at the Einstein radius. 



3.1 Mass profiles from gravitational lensing 

Strong and weak lensing data are used in combination 
to construct NFW profiles for the distant clusters A2390 
(Natarajan et al. 2005) and MS1358 (Hoekstra et al. 1998). 
For MS 1358 and A2390 the lensing mass profiles are found 
to be in excellent agreement with the X-ray mass profile 
(Arabadjis et al. 2002; Allen et al. 2001). In this work, we 
use the NFW fits to the total mass for all 5 clusters. 



3.2 Mass profiles from X-ray data 

We use NFW mass profiles derived from X-ray observations 
for both the nearby and distant sample clusters. Table 1 
presents the published fit parameters and error bars for all 
five clusters. 

On large scales, the X-ray profiles of most of our clus- 
ters - A2199, A496, A576 and MS1358 - appear smooth and 
circular, without any indications of recent merging activity 
(Markevitch et al. 1999, Kempner & David 2004, Arabadjis 
et al. 2002). However, high- resolution Chandra observations 
reveal substructure in the core regions of the nearby clus- 
ters. A496 has a cold front with a sharply decreasing surface 
brightness to the North (Dupke & White 2003). A576 has 
surface brightness edges probably due to gas stripped off a 
small merging subcluster in the centre (Kempner & David 
2004). The ICM of A2199 is interacting with lobes from the 
central radio source 3C 338 (Johnstone et al. 2002). The 
X-ray image of A2390 looks smooth but slightly elongated 
in different directions depending on radius. The central re- 
gion shows some substructure - a larger and a smaller X-ray 
peak, suggesting that the cluster is not fully relaxed after 
recent merger activity, although it appears to be in hydro- 
static equilibrium (Allen et al. 2001). 

The X-ray data provide estimates of the gas tempera- 
ture profiles needed to calculate the mass profiles and also 
f3x(r) (§6). Temperature profiles for the nearby clusters are 
available from both Chandra and ASCA/ROSAT observa- 
tions. The Chandra data have better spatial resolution than 
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Table 1. Parameters for the total mass profile from X-ray data, fitted to an NFW model with the reported 
errors in the fit parameters and the consequent errors on the total mass. The errors reported are l-cr errors. 



Cluster z r s , kpc c M(r = 0.5 Mpc), Mq Reference 



A2199 


0.0309 


129±i5 


fO 


(1.33 ± 0.15) x 10 14 


Markcvitch ct al. 1999 


A496 


0.032 


257±30 


6 


(1.66 ± 0.18) x 10 14 


Markevitch et al. 1999 


A576 


0.03829 


300±40 


4.5 


(1.11 ± 0.15) x 10 14 


Panes et al. 2000 


A2390 


0.2284 


6291™° 


3.3 


(3.26 ± 0.98) x 10 14 


Allen et al. 2001 


MSf358 


0.3289 


109t^ 6 


8.4 


(1.42 ± 0.28) x 10 14 


Arabadjis ct al. 2002 



the ASCA/ROSAT data, but extend to smaller radii (typ- 
ically about 0.2 Mpc as opposed to approximately 1 Mpc 
for the ASCA data). Chandra temperature data for A2199, 
A496 and A576 are obtained from Johnstone et al. (2002), 
Dupke & White (2003), and Kempner & David (2004), re- 
spectively. Johnstone et al. found the temperature profile 
of A2199 to be spherically symmetric to a good approxi- 
mation. Consequently, they binned the data into concentric 
annular rings, and fitted them using a single-temperature 
MEKAL model for each annulus (see Johnstone et al. 2002 
for details). We fit the Johnstone temperature profile with 
a logarithmic function using a least-squares algorithm. The 
best fit is: T = 1.657 lg(r) + 1.035, where T is in keV and r 
is in kpc. 

Due to significant temperature anisotropy, Dupke & 
White (2003) divided A496 into two regions (SHARP and 
SMOOTH; see their Fig. lb) and obtained separate tem- 
perature profiles for each region using a WABS MEKAL 
spectral model. We fit the temperature profiles with linear 
functions using a least-squares algorithm. The best fit to 
the SHARP data is T = 0.0475r + 1.75, where T and r 
are again in keV and kpc, respectively. The linear fit to the 
SMOOTH data, which is somewhat less perfect, produces 
very similar results for f3x- Therefore, we use the param- 
eters for the SHARP data hereafter. Finally, Kempner & 
David (2004) subdivided A576 into three sectors (see their 
Fig. 2). Their combined data are consistent with a constant 
temperature of 4 keV for the entire range of radii considered. 
Fig. 1 shows the temperature data and fits for A2199, A496 
and A576. All temperature profiles presented in this paper 
are three-dimensional (deprojected) profiles. Markevitch et 
al. (1998) showed that this model provides a good fit to an 
ensemble of 30 clusters observed with ASCA. De Grandi & 
Molendi (2002) found similar results for the average cluster 
temperature profile of 21 clusters observed with BeppoSAX. 
Because both of these satellites have poor angular resolution, 
the accuracy of these results has been controversial. How- 
ever, recent high-resolution Chandra measurements of tem- 
perature profiles in 11 clusters confirm that cluster temper- 
ature profiles do decline in the range 0.2-0.6 riso, although 
some clusters have "cooling cores" with smaller tempera- 
tures at radii <0.1riso (Vikhlinin et al. 2005). The apparent 
discrepancy between the ASCA and Chandra temperature 
profiles of A2199 and A496 in Figure 1 is thus most likely 
a result of the different spatial resolutions of the observa- 
tions. The ASCA temperature measurements are unreliable 
at very small radii due to poor spatial resolution, while the 
smaller field of view of Chandra prevents a direct comparison 
of the temperature profiles of these two clusters. The agree- 
ment between the average temperature profiles of 11 clus- 
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Figure 1. Temperature profiles for the nearby clusters. For A2199 
and A496, both the Chandra (solid line) and ASCA/ROSAT 
(dotted line) temperature fits are shown. The points are the ob- 
served data points. In the A576 plot the temperature data for 
two of the three cluster sectors (West and Southeast) described 
in Kempner & David 2004 arc shown along with our 4 keV fit 
line. 



ters observed with Chandra (Vikhlinin et al. 2005) and 16 
clusters observed with XMM-Newton (Piffaretti et al. 2005) 
with earlier ASCA and BeppoSAX results indicates that the 
ASCA temperature profiles are reliable at radii >0.1 Mpc. 
At smaller radii, the Chandra measurements are more reli- 
able, although the complex properties of the ICM may in- 
dicate that the ICM is not in equilibrium at these radii. 
The well-behaved declining temperature profiles at larger 
radii, however, are consistent with an ICM in approximate 
hydrostatic equilibrium. Because we wish to use the ICM 
to estimate the mass profile (and because the radial range 
more closely matches that of the cluster galaxies), we use 
the ASCA temperature profiles of A2199 and A496 for our 
dynamical analysis. 

The ASCA/ROSAT temperature data for A2199 and 
A496 are obtained from Markevitch et al. (1999). The au- 
thors fitted the data with polytropic models: 
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Figure 2. Temperature profiles from Chandra data for the dis- 
tant clusters. The points and solid lines indicate the observational 
data and fits, respectively. 
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where for A2199 r x = 95.7 kpc, a = 0.636, and 7 = 1.17, 
and for A496 r x = 178 kpc, a = 0.700, and 7 = 1.24. 
The ASCA data and fits for the nearby sample are shown in 
Fig. 1. For A576, Rines et al. (2000) find a temperature of 
3.77 ± 0.10 keV, and White (2000) finds 4.02 ± 0.07 keV, in 
good agreement with the Chandra value, suggesting that the 
cluster is approximately isothermal. The ASCA temperature 
profile of White (2000) is consistent with isothermal to a 
radius of 15' = 0.9Mpc. 

The temperature profiles of the distant clusters are 
based on Chandra data. For A2390, Allen, Schmidt & Fabian 
(2001) find that the temperature is well described by: 



T{r) 

72500 



To + Ti 



1 + (l/l c )l ' 



(12) 



where x = r/r 25 oo, T 25 oo = keV, r 25 oo = 0.64±g;J| 

Mpc, To = 0.40 ±0.02, Ti = 0.61 ±0.07, x c = 0.087 ±0.011, 
and r) = 1.9 ± 0.4. 

The temperature profile of MS1358 is consistent with an 
isothermal temperature of Tx = 7.16 ±0.10 keV (Arabadjis 
ct al. 2002). Fig. 2 shows the temperature data and fits for 
the distant sample. 



4 VELOCITY DISPERSION PROFILES 

The velocity data for our nearby clusters comes from the 
CAIRNS project, which combines new redshifts with those 
from many earlier studies. Optical and X-ray mass estimates 
of A576 disagree (Mohr et al. 1996, Rines et al. 2000), so it 
is of great interest to determine if physically reasonable or- 
bital distributions can reconcile this conflict. Markevitch et 
al. (1999) proposed that A496 and A2199 are good exam- 
ples of relaxed clusters based on their circularly symmetric 
X-ray contours and their well-behaved temperature profiles. 
Despite its quiescent central region, A2199 is surrounded by 
several infalling X-ray groups and has an asymmetric galaxy 
distribution (Rines et al. 2001, 2002). 

We calculate velocity dispersion profiles (VDPs) for 
these three clusters using galaxies defined to be members 
from the caustic technique. In redshift space, cluster infall 
regions exhibit dense envelopes in redshift-radius diagrams. 



These envelopes, known as caustics, provide a straight- 
forward membership classification. The use of caustics for 
denning cluster membership is similar to the "shifting gap- 
per" technique (Fadda et al. 1996), but is less sensitive to 
the density of the galaxy sample. We use the caustics found 
by Rines et al. (2003) to define membership for the nearby 
clusters. VDPs calculated with only galaxies brighter than 
M K = -23.5 (Rines et al. 2004) in the Two-Micron All-Sky 
Survey (Nikolaev et al. 2000) are consistent with those used 
here, suggesting that bright and faint galaxies have similar 
orbital distributions and velocity biases. We calculate the 
velocity dispersion in bins of 25 galaxies; changing the bin 
size does not significantly affect the VDPs. The typical red- 
shift uncertainties are ~30 km s~\ so these uncertainties 
contribute only a small amount of uncertainty to the VDPs. 
We do not attempt to subdivide the galaxy sample to obtain 
different orbital solutions for different sub-samples (Biviano 
& Katgert 2004). 

Like the total mass, the galaxy number density of the 
nearby clusters is modeled with an NFW profile. Rines et 
al. (2004) calculate the number density profiles of galaxies 
in A2199, A496 and A576 from complete, Jf 3 -selected spec- 
troscopic samples. The galaxy samples are complete to 1.5-2 
magnitudes fainter than M* and membership is defined from 
spectroscopic redshifts using the caustic technique (Diaferio 
1999). Rines et al. (2004) fit the observed number density 
profiles to NFW models and find scale radii of r Stgxy =HA 
kpc, 500 kpc, and 1.07 Mpc for A496, A576, and A2199 
respectively. Because these profiles are derived from com- 
plete, /^-selected spectroscopic samples, there is very little 
uncertainty in the profiles due to incompleteness. The over- 
all normalization of the galaxy density (and thus the value 
of c) is immaterial since it cancels out in the calculation of 
/?orb(r) (see Section 2). The total mass, line-of-sight velocity 
dispersion and galaxy number density profiles are all the in- 
gredients needed to calculate f3 OTh (r), as explained in Section 
2. 

For the two distant clusters, we again use the caustics 
to define cluster membership (Diaferio et al. 2005). Red- 
shift data are collected from the literature, with the major- 
ity coming from Yee et al. (1996) [A2390] and Fabricant 
ct al. (1991), Yee et al. (1998), and Fisher et al. (1998) 
[MS 1358]. 

The VDPs for all five clusters are fitted with polynomial 
functions using a least-squares algorithm. The lowest degree 
polynomial producing a good fit (typically \ 2 /DOF < 1) is 
used. The polynomial fits are applied in the radial range < 
r < r m ax, where r max is the maximum radius for which VDP 
data are available. For r max < r < 5 Mpc, ai os is set equal 
to the value at r max , and beyond that it is set equal to zero. 
We exclude the innermost data points for A496 and A2390 
because including them requires higher-order polynomials 
to yield acceptable fits to the VDP data; these polynomials 
lead to unphysical solutions for f3 orb (r). Figure 3 shows the 
VDP data and fits for the five clusters. 

The galaxy number densities of A2390 and MS 1358 are 
well-described by a Hernquist profile: 



v g {r) oc 



1 



(r/r c )(l+r/r c )3 



(13) 



with r c = 1.42 Mpc and r c = 1.08 Mpc, respectively (Carl- 
berg et al. 1997). 
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Figure 3. Velocity dispersion profiles and polynomial fits for all 
five clusters. 



The CNOC1 group estimate the galaxy number den- 
sity profiles of A2390 and MS 1358 from their large spec- 
troscopic samples including a completeness correction for 
galaxies without redshifts which depends on both position 
and magnitude (Carlberg et al. 1997). They find that the 
number density profiles are well fit by a Hernquist profile: 
[eqn. 13] with r s , gxy =lA2 and 1.08 Mpc for A2390 and 
MS1358 respectively. 



5 RESULTS FOR THE VELOCITY 
ANISOTROPY PARAMETER 

The results for /3 or b for all five clusters are shown in the left- 
hand panels of Fig. 4. We find diversity in orbital parameters 
in our sample. Two of the clusters A2199 and A496 in the 
nearby sample appear to have predominantly tangential or- 
bits in the inner region (out to 1 Mpc) as indicated by the 
behaviour of (3 OT b(r), whereas A576 appears to have mostly 
radial orbits. The distant clusters seem to be dominated by 
radial orbits. Computing fix using the observed tempera- 
ture profile and the cr^(r) computed from the anisotropic 
Jeans equation, we find that MS1358 and A576 appear to 
deviate from hydro-static equilibrium (fix 7^ 1) from the 
inner regions to all the way out, and the rest seem to be 
consistent with being by-and large in hydro-static equilib- 
rium. For A576 and MS1358, fix > 1 indicating less energy 
in the gas compared to the galaxies, whereas for A2199, 
A496 and A2390 fix < 1. For A2199 and A496, the (3 X pro- 
files computed from the ASCA/ROSAT temperature data 
are slightly higher than those from the Chandra data and 
have a positive slope. Henceforth, we focus on the former 
because of the larger spatial extent necessary for compar- 
ison with the numerical simulations. We show in Figure 5 



Figure 4. Velocity anisotropy parameter /3 or b and energy ra- 
tio fix- The top- left panel shows /3 or b for the nearby clusters: 
A2199 (red curve), A496 (blue curve) and A576 (black curve). 
The top-right panel shows fix f° r the nearby clusters computed 
using ASCA/ROSAT and Chandra temperature data. Colors are 
as for the left panel: solid red - A2199 {ASCA/ROSAT), solid blue 
- A496 (ASCA/ROSAT), black - A576 (ASCA and Chandra), red 
dashed - A2199 (Chandra) and blue dashed - A496 (Chandra). 
The bottom-left panel shows /3 or b for the high-redshift clusters. 
Lines are as follows: green curve - MS1358 and cyan curve - A2390. 
Finally, the bottom-right panel shows fix for the distant clusters 
(lines are as in the bottom-left panel). The vertical ticks show 
the maximum extent of the X-ray temperature data for the left 
panels. For the right panel the vertical ticks show the extent to 
which the mass profile is fitted. The thick line marks /3 or b = 
and fix = 1. 

the effect of changing the normalization of the total VDPs 
by ±10%. 



6 COMPARISON WITH COSMOLOGICAL 
SIMULATIONS 

We now compare the observational results for /3 or b, fix 
and S(r) with those of high-resolution cosmological sim- 
ulations based on the "concordance" flat ACDM model: 
O m = 1 - fi A = 0.3, Q b = 0.021 h~ 2 , h = 0.7 and cr 8 = 0.9, 
where the present-epoch Hubble constant is defined as 100 h 
km s -1 Mpc -1 , and as is the power spectrum normalization 
on 8 h~ J Mpc scale. The simulations are performed using 
the Adaptive Refinement Tree (ART) N-body+gasdynamics 
code (Kravtsov, Klypin & Khokhlov 1997; Kravtsov 1999; 
Kravtsov, Klypin & Hoffman 2002), an Eulerian code with 
adaptive refinement in space and time and non-adaptive re- 
finement in mass. 

To set up initial conditions we first run a low-resolution 
simulation of 80 h~ 1 Mpc and 120 /i -1 Mpc boxes and select 
six clusters with virial masses ranging from r;6x 10 14 to 
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1.6 x 10 15 h' 1 Mg. The largest one has a virial mass ol 
1.6 x 10 15 h~ x Mq, the second massive cluster has a mass of 
1.1 x 10 lj hT 1 Mq, and the other four clusters have masses 
of ~ 7 — 10 x 10 14 Mq. The corresponding virial radii, 
defined as radii enclosing the overdensity of 337 with respect 
to the mean density of the Universe, for the six clusters are: 
2.36, 2.10, 2.01, 1.93, 1.92, and 1.78 hT 1 Mpc. 

The initial conditions are set using multiple-mass par- 
ticle technique retaining the previous large-scale waves in- 
tact but including additional small-scale waves, as described 
by Klypin et al. (2001). The re-sampled lagrangian region 
of each cluster, corresponding to a sphere of (1.5 — 5)7? v ir 
around it at Z = 0, is then simulated with high dynamic 
range. The actual spatial resolution of the simulations is ~ 5 
h' 1 kpc. The mass resolution (i.e., the dark matter particle 
mass) is m p ~ 9.1 x 10 s hT 1 Mq for the all but the most mas- 
sive cluster and 2.7 x 10 s hT 1 Mq for the most massive clus- 
ter. The clusters thus have million or more dark matter par- 
ticles within their virial radii. The simulations follow dissi- 
pationless dynamics of dark matter particles and gasdynam- 
ics of the baryonic component and include several processes 
crucial to galaxy formation: star formation, metal enrich- 
ment and thermal feedback due to supernovae types II and 
la, self-consistent advection of metals, metallicity-dependent 
radiative cooling and UV heating due to cosmological ion- 
izing background. Stellar feedback on the surrounding gas 
includes injection of energy and metals via stellar winds and 
supernovae as well as secular mass loss (see Kravtsov, Nagai 
& Vikhlinin 2005, for details). 

For each cluster simulation, we identify the main clus- 
ter and its galaxies using a variant of the Bound Density 
Maxima (BDM) algorithm using dark matter particles only. 
The details of the algorithm and parameters used in the 
halo finder can be found in Kravtsov et al. (2004). The 
main steps of the algorithm are identification of local density 
peaks (potential halo centers) and analysis of the density dis- 
tribution and velocities of the surrounding particles to test 
whether a given peak corresponds to a gravitationally bound 
clump. More specifically, we construct density, circular ve- 
locity, and velocity dispersion profiles around each center 
and iteratively remove unbound particles using the proce- 
dure outlined in Klypin et al. (1999). We then construct fi- 
nal profiles using only bound particles and use these profiles 
to calculate proper ties of halos, such as the circular velocity 
profile Vcirc(r') = \l GM{< r)/r and compute the maximum 
circular velocity Vn lax . For halos located within the virial 
radius of a larger host halo {the subhalos), we define the 
outer boundary at the truncation radius, r t , at which the 
logarithmic slope of the density profile constructed from the 
bound particles becomes larger than —0.5 as we do not ex- 
pect the density profile of the CDM halos to be flatter than 
this slope. For each system we estimate the stellar mass, 
M„, gas mass and total mass (dark matter, stars, and gas) 
within the truncation radius. 

Once the galaxies are identified, the radial and tan- 
gential components of the velocity dispersion, ay and at, 
of the dark matter, gas and galaxies are measured in ra- 
dial bins centered on the cluster potential minimum after 
subtracting the peculiar velocity of the cluster, defined as 
the mass-weighted bulk velocity of dark matter within the 
cluster core. Only galaxies with masses M, > 1 x 10 9 Mq 
(or Vmax > 80 kms -1 ) are used in this calculation. This 
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Figure 5. Velocity anisotropy /3 or b - comparison between obser- 
vations and numerical simulations. The error bars on each curve 
shows the error in /3 or b due to changing <r; os by ±10%, evalu- 
ated at the given radius. The points represent the six simulated 
clusters. 



value for the threshold provides a reasonably large sample 
of galaxies, while not compromising the numerical resolu- 
tion of the results. The velocity anisotropy parameter fJ OI b 
is computed from the components of the velocity dispersion 
using the definition in Section 2. In addition, fix is com- 
puted using the galaxy of and the X-ray temperature of the 
gas. The measured temperatures are gas mass weighted is 
calculated assuming Chandra energy response in the 0.5-7 
keV band. The temperature profiles in the simulations are 
extracted from the 3D distribution and have negligible mea- 
surement errors. 

Unfortunately, the numerical results are rather noisy 
in the central cluster regions due to the smaller number of 
galaxies there. To allow a better comparison with the numer- 
ical simulations, we extend the observed /3 rb and fix profiles 
to 2 Mpc by increasing the integration limit Tit from 3.5 to 
7 Mpc. Beyond 1.5 Mpc, there is a small amount of scat- 
ter in the (3 OI b profile due to numerical effects. Increasing 
Rt beyond 7 Mpc slows down the calculation considerably 
without reducing this scatter much further, so 7 Mpc is es- 
tablished to be the optimal value for the large R limit for 
the purposes of this calculation. 

The detailed results of the comparison with the simula- 
tions are discussed in the section below. 



7 DISCUSSION AND CONCLUSIONS 

The statistical uncertainties in the individual data points for 
the VDPs and the number density profiles are ~20%; the 
uncertainties in the mass profiles are ~10%. To illustrate 
the impact these uncertainties have on the derived galaxy 
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Figure 6. Energy ratio fix - comparison between observations 
and numerical simulations. The /3x's are computed from the /3 or b 
using ASCA/ROSAT temperature data for the nearby clusters 
and Chandra data for the distant ones. The error bars on each line 
show the error in fix d ue to changing cr ios by ±10%, evaluated at 
the given radius. The points represent the six simulated clusters. 

orbits, we recalculate the orbital distribution j3 rb(r) after 
renormalizing the entire VDP by ±10% (the uncertainty in 
the total VDP is smaller than the individual bins). Because 
the /3 or b(r) profiles are much more sensitive to changes in 
VDP than either v(r) or M(r), these uncertainties are rep- 
resentative of the total uncertainty from the observations 
(numerical errors are significantly smaller than the observa- 
tional uncertainties). Figures 5 and 6 show the magnitude 
of these uncertainties; they are comparable to the scatter 
seen in the simulated clusters. The derived orbital profiles 
agree well with the simulated clusters, although the uncer- 
tainties in both observations and simulations are still large 
enough to make it difficult to distinguish mildly radial orbits 
from mildly circular orbits. As the observations and simu- 
lations improve, more detailed treatments of the systematic 
uncertainties will be required. The current study demon- 
strates the feasibility of this technique and the consistency 
of various mass estimators. That is, we successfully derive 
physical orbit models for the cluster galaxies assuming that 
they trace the mass profile inferred from independent X-ray 
or lensing data. We further show that these orbital profiles 
closely resemble those of cluster galaxies in simulations 

In our cluster sample, we find a diversity of galaxy or- 
bits. All nearby clusters appear to have primarily tangential 
orbits. The distant cluster MS1358 exhibits radial orbits 
while A2390 has predominantly tangential orbits at small 
radii but radial orbits at large radii (r > 1 Mpc). The sta- 
tistical uncertainties in the observed quantities can lead to 
large uncertainties in f) OI b' the statistical ~ la uncertain- 
ties are ~ 20% for the VDPs, ~ 10% for the mass profiles, 
and « 20% for the galaxy number densities. As an illustra- 



tion, the errors plotted in Fig. 6 and Fig. 7 are calculated 
by changing ai os by ±10%, while keeping the total mass and 
galaxy number density fixed. Because /3 or b is much more sen- 
sitive to changes in the VDP than to the other two quanti- 
ties, these errors are representative of the overall uncertainty 
in f3 OT b and fix ■ Numerical errors are generally within a few 
per cent and much smaller than the observational uncertain- 
ties. At larger radii, where observational constraints (lensing 
and X-ray data) need to be extrapolated, /3 or b is not directly 
constrained by the data. 

We find departures from hydrostatic equilibrium in all 
of our clusters except A2390 (especially when the X-ray mass 
is used to compute /3 rb)- However, these deviations are gen- 
erally small (0.6 < fix < 1-6 excluding the CTiandra-based 
results for A2199 and A496) and comparable to the mar- 
gins of error (see Fig. 6). This implies that the cluster mass 
estimates from X-ray observations, which are based on the 
assumption of hydrostatic equilibrium, may still be valid. 
A2199, A496 and A2390 have fix < 1, indicating more en- 
ergy in the ICM than in the galaxies, whereas A576 and 
MS1358 have fix > 1, indicating less energy in the gas. 

There is general agreement in shape, magnitude and 
range between the observed and simulated /3 or b profiles (Fig. 
5), especially for r > 0.5 Mpc, where the simulations are re- 
liable due to sufficient galaxy statistics. This agreement im- 
plies that the anisotropic Jeans equation with its assump- 
tions of spherical symmetry is applicable to observational 
data. The simulated fix profiles also have similar shapes 
and ranges to the observed ones, but their values are ap- 
proximately a factor of two larger (Fig. 6). The reason for 
this discrepancy is not clear. Our results are consistent with 
theoretical expectations that clusters are dynamically com- 
plex as they are still in the process of assembling. The dif- 
ference between the inferred orbits is curious because it sug- 
gests that the moderate redshift clusters are dynamically 
younger than the nearby clusters. Due to the small sample 
size, however, it would be premature to conclude that clus- 
ter galaxy orbits evolve dramatically between 2=0.2-0.3 and 
the present. If this effect is real, it could result from the ob- 
served evolution in the infall rate of galaxies (Ellingson et 
al. 2001). 

There are several directions in which the present work 
could be extended. First, one could use a larger cluster sam- 
ple to verify the statistical significance of the results pre- 
sented in this paper - e.g. whether the differences in orbital 
structure between the nearby and distant clusters are a co- 
incidence due to our small samples or a more general fact. 
In addition, better observational data, especially more pre- 
cise cluster VDPs, would constrain the velocity anisotropy 
much further. The numerical simulations could also be re- 
fined in order to resolve better the central cluster regions. 
New sources of entropy injection could be explored with sim- 
ulations that include processes such as heating by quasars 
and AGNs, thermal conduction and turbulent mixing. 
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